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Estimate of 13{B X^-f) at 0{ct'^) 
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Combining our results for various 0(0^) corrections to the weak radiative _B-meson decay, we 
are able to present the first estimate of the branching ratio at the next-to-next-to-leading order in 
QCD. We find B{B X^j) = (3.15 ± 0.23) x 10~* for > 1.6 GeV in the B-meson rest frame. 
The four types of uncertainties: nonperturbative (5%), parametric (3%), higher-order (3%) and 
jTic-interpolation ambiguity (3%) have been added in quadrature to obtain the total error. 

PACS numbers: 12.38.Bx, 13.20.He 



The inclusive radiative _B-meson decay provides im- 
portant constraints on the minimal supersymmetric stan- 
dard model and many other theories of new physics at the 
electroweak scale. The power of such constraints depends 
on the accuracy of both the experiments and the stan- 
dard model (SM) calculations. The latest measurements 
by Belle and BABAR are reported in Refs. [E i]. The 
world average performed by the Heavy Flavor Averaging 
Group [3] for > 1.6 GeV reads 




(3.55 ± 0.24 tofo ± 0.03) x 10"^ 



(1) 



The combined error in the above result is of the same 
size as the expected 0{al) next-to-next-to-leading or- 
der (NNLO) QCD corrections to the perturbative de- 
cay width r(6 Ar|"'''*°"7), and larger than the known 
nonperturbative corrections to the relation T{B 
Xs-f) ~ r(6 ^ XP'^'-*°"7) Thus, calculating the 

SM prediction for the 6-quark decay rate at the NNLO is 
necessary for taking full advantage of the measurements. 

Evaluating the ©(a^) corrections to B{b Xp^'^"""-/) 
is a very involved task because hundreds of three-loop 
on-shell and thousands of four-loop tadpole Feynman di- 
agrams need to be computed. In a series of papers [3]- 
[14| . we have presented partial contributions to this en- 
terprise. The purpose of the present Letter is to combine 
all the existing results and obtain the first estimate of 
the branching ratio at the NNLO. We call it an estimate 
rather than a prediction because some of the numeri- 
cally important contributions have been found using an 
interpolation in the charm quark mass, which introduces 
uncertainties that are difficult to quantify. 



FIG. 1: Sample LO diagram for the 6 ^ 57 transition. 



Let us begin with recalling that the leading-order (LO) 
contribution to the considered decay originates from one- 
loop diagrams in the SM. An example of such a diagram 
is shown in Fig. [T] Dressing this diagram with one or 
two virtual gluons gives examples of diagrams that one 
encounters at the next-to-leading order (NLO) and the 
NNLO. In addition, one should include diagrams describ- 
ing the bremsstrahlung of gluons and light quarks. 

An additional difficulty in the analysis of the con- 
sidered decay is the presence of large logarithms 
(a^ In M^/m^)" that should be resummed at each or- 
der of the perturbation series in as- To do so, one em- 
ploys a low-energy effective theory that arises after de- 
coupling the top quark and the heavy electroweak bosons. 
Weak interaction vertices (operators) in this theory are 
either of dipole type {saf'^hF^,^, sa^"'T''hGl,,) or con- 
tain four quarks ([sF^] [gr'g]). Coupling constants at 
these vertices (Wilson coefficients) are first evaluated 
at the electroweak renormalization scale /ig ~ , Mw 
by solving the so-called matching conditions. Next, 
they are evolved down to the low-energy scale ~ frib 
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according to the effective theory renormahzation group 
equations (RGE). The RGE are governed by the oper- 
ator mixing under renormahzation. Finally, one com- 
putes the matrix elements of the operators, which in 
our case amounts to calculating on-shell diagrams with 
single insertions of the effective theory vertices. 

A summary of the B — > Xa^ calculation status be- 
fore the beginning of our project can be found, e.g., in 
Ref. [Hi. At the NNLO level, the dipole and the four- 
quark operators need to be matched up to three and two 
loops, respectively. Renormahzation constants up to four 
loops must be found for 6 — > 57 and b sg diagrams with 
four-quark operator insertions, while three-loop mixing is 
sufficient in the remaining cases. Two-loop matrix ele- 
ments of the dipole operators and three-loop matrix ele- 
ments of the four-quark operators must be evaluated in 
the last step. 

Three-loop dipole operator matching was found in 
Ref. [1]. The necessary three-loop mixing was calcu- 
lated in Ref. @. The four-loop mixing was evaluated in 
Ref. [1^ . Two- loop matrix element of the photonic dipole 
operator together with the corresponding bremsstrahlung 
was found in Refs. [HI and recently confirmed in 
Ref. Three- loop matrix elements of the four-quark 

operators were found in Ref. within the so-called 
large-/3o approximation. A calculation that goes be- 
yond this approximation by employing an interpolation 
in the charm quark mass rUc has just been completed in 
Ref. 

With all these results at hand, we are ready to present 
the first estimate of the B ^ Xg^i branching ratio at 
0{al). It reads [H 



B{B Xa-f) = (3.15 ± 0.23) x 10"'', 



(2) 



for Ej > 1.6 GeV in the i?-meson rest frame. The four 
types of uncertainties: nonperturbative (5%), parametric 
(3%), higher-order (3%) and mc-interpolation ambiguity 
(3%) have been added in quadrature in Eq. 

The central value in Eq. ^ was obtained for /io ~ 
160 GeV, fit = 2.5 GeV and fic = 1.5 GeV. The latter 
quantity stands for the charm mass MS renormahzation 
scale that is allowed to be different from /i^. The branch- 
ing ratio dependence on each of the three scales is shown 
in Fig. [21 Once one of them is varied, the remaining 
two are fixed at the values that have been mentioned 
above. The reduction of the renormahzation scale de- 
pendence at the NNLO is clearly seen. The most pro- 
nounced effect occurs for Hc that was the main source of 
uncertainty at the NLO. (The LO results are nic- and 
thus /ic-independent.) The current uncertainty of ±3% 
due to higher-order [©(a^)] effects is estimated from the 
NNLO curves in Fig. d 

The reference value oi fib — 2.5 GeV that we have 
chosen is roughly twice smaller than in the previous LO 
and NLO analyses. Given the stability of the NNLO 
result for large values of /Zb, we do not underestimate 
any uncertainty from that region. Furthermore, be- 
cause the center-of-mass energy ttib — 5.3 GeV gets 
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FIG. 2: Renormalization scale dependence of B(B — > Xs'y) in 
units 10"* at the LO (dotted lines), NLO (dashed lines) and 
NNLO (solid lines). The plots describe subsequently the depen- 
dence on the matching scale fio, the low-energy scale fit, and 
the charm mass renormalization scale /ic. 



distributed among various partons, the reference value 
of /ift = 2.5 GeV seems reasonable. Lower values of fit 
have an advantage of making /ic-stabilization more ef- 
ficient because the NNLO logarithm that compensates 
/ic-dependence of the NLO amplitude comes multiplied 
by asinb)- 

The ±3% uncertainty that is assigned to the nic- 
interpolation ambiguity has been estimated studying by 
how much the NNLO branching ratio depends on vari- 
ous interpolation assumptions. More details on this point 
and other elements of the phenomenological analysis (in- 
cluding the input parameters) can be found in Ref. 14]. 

As far as the parametric uncertainties are concerned, 
the dominant ones come from as{Mz) (±2.0%) and the 
measured semilcptonic branching ratio B{B — > XcBD) 
(±1.6%) to which we normalize. The third-to-largest 
uncertainty (±1.1%) is due to the correlated errors in 
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nT-cifnc) and the semileptonic phase-space factor 



C = 



Vub 



T[B ^ X^ev] 



(3) 



The factor C has been determined in Ref. together 
with nic^mc) from a global fit to the semileptonic data. 
If the normalization to B{B — > XcsD) was not applied in 
the B — > Xs"f calculation, the error due to mc{mc) would 
amount to ±2.8%. At the same time, one would need to 
take into account uncertainties in and the Cabibbo- 
Kobayashi-Maskawa factor |V(*Vt6p, each of which ex- 
ceeds ±3%. 

The nonperturbative uncertainty in Eq. ([2]) is due to 
matrix elements of the four-quark operators in the pres- 
ence of one gluon that is not soft (Q^ ~ m^,TOfcA, where 
A ~ Aqcd)- Unknown nonperturbative corrections to 
them scale like ofsA/mh in the limit rric ^ Tib/Z and like 
UsK^ /m^ in the limit nic ^ mb/2. Because rric < mb/2 
in reality, asA^/mi, should be considered as the quan- 
tity that sets the size of such effects. Consequently, a 
±5% nonperturbative uncertainty has been assigned to 
the result in Eq. Q . This is the dominant uncertainty at 
present. Thus, a detailed analysis of such effects would 
be more than welcome. So far, no published results on 
this issue exist. Even lacking a trustworthy method for 
calculating such effects, it might be possible to put rough 
upper bounds on them that could supersede the current 
guess-estimate of ±5%. Nonperturbative corrections to 
inclusive B — > Xd.s^ decays that scale like A/rub may 
arise when the 6-quark annihilation vertex does not co- 
incide with the hard photon emission vertex; see, e.g., 
Ref. or comments on B ^ Xdj in Sec. 2 of Ref. [s]. 

The NNLO central value in Eq. ^ differs from some 
of the previous NLO predictions by between 1 and 2 error 
bars of the NLO results. Because those error bars were 
obtained by adding various theoretical uncertainties in 
quadrature, such a shift is not improbable, similarly to 
shifts by less than 2a in experimental results. The shift 
from the NLO to the NNLO level diminishes with low- 
ering the value of fic, which has motivated us to use the 
relatively low fic = 1.5 GeV as a reference value here. 

The NNLO results turn out to be only marginally de- 
pendent on whether one follows (or not) the approach 
of Ref. where the top-quark contribution to the de- 
cay amplitude was calculated separately and rescaled by 
quark mass ratios to improve convergence of the pertur- 
bation series. Although the top contribution alone in- 
deed behaves better also at the NNLO level when such 
an approach is used, the charm quark contribution (to 
which no rescaling has been applied in Ref. [17]) does 
not turn out to be particularly stable beyond the NLO. 
Consequently, in the derivation of Eq. ^ and Fig. [21 wc 
have used the simpler method of treating charm and top 
sectors together. 

Our result in Eq. ([2]) has been obtained under the as- 
sumption that the photonic dipole operator contribution 
to the integrated spectrum below 1.6 GeV is well ap- 
proximated by a fixed-order perturbative calculation (see 
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FIG. 3: B{B Xs"f) as a function of the charged Higgs boson 
mass in the THDIVI II for tan (3 = 2 (solid lines). The dashed and 
dotted lines show the SM and experimental results, respectively 
(see the text). 



Note added). For lower values of the photon energy cut, 
the following numerical fit can be used: 



/ B{E., > Eq) \ 
\B{E^ > 1.6 GeV) / 



l + 0.15x-0.14a;^ (4) 



fixed 
order 



where x = 1 — i?o/(l-6 GeV). This formula coincides with 
our NNLO results up to ±0.1% for Eq S [1.0, 1.6] GeV. 
The error is practically i?o-independent in this range. 

In the remainder of this Letter, we shall update the 
B ^ Xsj constraints on the charged Higgs boson mass 
in the two-Higgs-doublet-model II (THDM II)J18]. The 
solid lines in Fig.[3]show the dependence of B{B Xg'y) 
on this mass when the ratio of the two vacuum expecta- 
tion values, tan/3, is equal to 2. The dashed and dotted 
lines show the SM (NNLO) and the experimental results, 
respectively. In each case, the middle line is the cen- 
tral value, while the other two lines indicate uncertainties 
that one obtains by adding all the errors in quadrature. 

In our THDM calculation, matching of the Wilson co- 
efficients at the electroweak scale is complete up to the 
NLO 19], but the NNLO terms contain only the SM con- 
tributions (the THDM ones remain unknown) . In conse- 
quence, the higher-order uncertainty becomes somewhat 
larger. This effect is estimated by varying the matching 
scale fj,o from half to twice its central value. It does not 
exceed ±1% for the Mh+ range in Fig.[3l 

Even though the experimental result is above the SM 
one, the lower bound on Mh+ for a generic value of 
tan [3 remains stronger than what one can derive from 
any other currently available measurement. If all the 
uncertainties are treated as Gaussian and combined in 
quadrature, the 95% (99%) CL bound amounts to around 
295 (230) GeV. It is found for tan (3 oo but stays prac- 
tically constant down to tan/3 ~ 2. For smaller tan/3, 
the branching ratio and the bound on Af^+ increase. 

The contour plot in Fig. [3] shows the dependence of 
the Mjj+ bound on the experimental central value and 
error. The current experimental result ([1]) is indicated by 
the black square. Consequences of the future upgrades in 
the measurements will easily be read out from the plot, 
so long as no progress on the theoretical side is made. Of 
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FIG. 4: The 95% CL lower bound on as a function of the 

experimental central value (horizontal axis) and error (vertical 
axis). The experimental result from Eq. dTJ is indicated by the 
black square. The contour lines represent values that lead to the 
same bound. 

course, the derived bounds should be considered illustra- 
tive only because they depend very much on the theory 
uncertainties that have no statistical interpretation. 

To conclude, we have provided the first estimate of 
B{B Xs-f) at 0{al). The inclusion of the NNLO 
QCD corrections leads to a significant suppression of the 
branching ratio renormalization scale dependence that 
has been the main source of uncertainty at the NLO. 



The central value is shifted downward with respect to 
all the previously published NLO results. It is now 
about Icr lower than the experimental average ([1]). The 
dominant theoretical uncertainty is currently due to the 
unknown O (asA/mb) nonperturbative effects. In the 
two-Higgs-doublet model II, the experimental results fa- 
vor a charged Higgs boson mass of around 650 GeV. 
The 95% C.L. bound for this mass amounts to around 
295 GeV if all the uncertainties are treated as Gaussian. 
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Note added — Recently, our results from Eqs. and 
(|4|) were combined in Ref. [l^l with perturbative cutoff- 
related corrections that go beyond a fixed-order calcula- 
tion 0,[2l|. Because these corrections ioi Eq < 1.6 GeV 
do not exceed our higher-order uncertainty of ±3%, we 
postpone their consideration to a future upgrade of the 
phenomenological analysis, where other contributions of 
potentially the same size are going to be included, too 
(see Sec. 1 of Ref. f22ll). 
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